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Abstract 

We propose a nonadiabatic approach to quantum anneahng, in which we repeat quantum an- 
neahng in nonadiabatic time scales, and collect the final states of many realizations to find the 
ground state among them. In this way, we replace the diffculty of long annealing time in adiabatic 
quantum annealing by another problem of the number of nonsidabatic (short-time) trials. The 
one-dimensional transverse-field Ising model is used to test this idea, and it is shown that nona- 
diabatic quantum annealing has the same computational complexity to find the ground state as 
the conventional adiabatic annealing does. This result implies that the nonadiabatic method may 
be used to replace adiabatic annealing to avoid the effects of external disturbances, to which the 
adiabatic method is more prone than the nonadiabatic counterpart. 
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I. INTRODUCTION 



Quantum annealing (QA) is an algorithm to solve optimization problems by the method 



of quantum mechanics and statistical physics 



-|8|. A combinatorial optimization problem 



can often be expressed as a problem to find the ground state of a spin system, where the cost 
function is described as the Hamiltonian ^H|. Most current research activities on QA are 



12 1 . The adiabatic theorem states 



based on the adiabatic theorem of quantum mechanics 
that, if the time evolution of the Hamiltonian is sufficiently slow, the quantum system stays 
in the instantaneous ground state when we start from the ground state. Therefore, if we 
choose the initial Hamiltonian to have a trivial ground state and the final Hamiltonian as the 
target Hamiltonian whose ground state we are to find, the solution is automatically reached 
after the adiabatic time evolution following the time-dependent Schrodinger equation. 

To describe the problem more quantitatively, let us consider the following time-dependent 
Hamiltonian 

Hit)=(^l-^-^Ho+^-V. (1) 

Initially (t = 0), the Hamiltonian H{0) equals to V with a trivial ground state, and at the 
final time t = t the Hamiltonian H{t) reduces to the target Hamiltonian Hq. t determines 
the time scale of a single realization of annealing. In order to apply the adiabatic theorem, 
the following inequality has to be satisfied: 

. maxo<,<.[(l(t)|^|0(t))] 

mino<i<.[A(t)2] ' ^ ^ 

where |0(t)) and express the instantaneous ground and first excited states of the 

Hamiltonian at time t, respectively, and A describes the energy gap between those states. 

In practice, the minimum energy gap mino<t<r[A(t)] often increases rapidly as the system 
size N increases, and correspondingly, the annealing time r diverges. Especially, for problems 
called NP hard 13| , the minimum energy gap decreases exponentially with N. This difficulty 
of adiabatic quantum annealing is our motivation to follow a nonadiabatic approach to 
quantum annealing. 

In the present paper, we propose a nonadiabatic method to find the ground state of the 
target Hamiltonian. For this purpose, we repeat the realizations of quantum annealing, and 
a single realization is done in a nonadiabatic time scale. We collect the final states of many 
realizations with the expectation that the ground state is found among them. In this way. 



we replace the difficulty of long annealing time in adiabatic quantum annealing by another 
problem of the number of nonadidabatic (short-time) trials. 

We show that nonadiabatic quatum annealing can succeed in the same time scale (com- 
putational complexity) as the adiabatic one for the one-dimensional transverse-field Ising 
model. Nonadiabatic QA has the advantage of an easier practical implementation because 
we do not have to tune the time-dependent parameters carefully to keep the system in the 
ground state. This fact also implies that the nonadiabatic method may be more robust 



against external distur 
nealing, Johnson et al 



nances. In fact, in an experimental implementation of quantum an- 



14| repeated annealing processes and measured the probability that 



the final system is in the ground state (See also 15|). 

This paper is organized as follows. In the next section, we formulate the idea explicitly 
for the transverse-field Ising chain and analyze the system analytically and numerically. The 
final section is for summary and discussion. 

In the whole of this paper, we use the units = h = 1. 



II. NONADIABATIC QUANTUM ANNEALING 

Let us introduce the one-dimensional transverse-field Ising model: 

f ^ / f\ ^ 

^w = --E-^^-m- 1-- E-^ (3) 

i=l ^ ^ i=l 

The Pauli matrix af'^ satisfies the periodic boundary condition c^^i = o"^'^. The first 
coupling term corresponds to the target Hamiltonian Hq in Eq. ([1]), and the second term is 
the initial Hamiltonian V . 



Following the standard procedure 16|,ll7|], the Hamiltonian is decomposed into the sum of 
the two-state space of each mode q = ivr/A^, iSvr/A^, . . . , ±(A^— l)7r/A^. The time-dependent 
Schrodinger equation is reduced to a simple equation 

d fUq\ f {1 - t/r) - (t/r) cos q (t/r)sing \ f'^i 



2^ M = 2 r ' ' - ^ ^ - ^ - (4) 

(^t\Vq/ \ (t/r)sing -{1 - t/r) + (t/r) cosql \Vg/ 

The adiabatic condition, applied to the smallest gap among various q values, turns out to 
be 



3 



The energy gap A,(t) takes ,ts mimmum a» = r/2 for a,ry 

According to Vitanov and Garraway [18], the dynamics of such a two- level quantum 
system is exactly solvable and the solution is expressed as a linear combination of parabolic 
cylinder functions Dy{z). However, we follow an approximate method, which is expected to 
be asymptotically correct in the limit ^ 1 and is easier to see explicitly. 

First, we show the energy spectra ie^ for q = O.Tvr, O.Svr, O.OSvr in Fig. [H It is observed 
that a state transition is likely to happen only around t = r/2, where the gap is very small 
for small q. This fact allows us to replace t/r in the off-diagonal elements in Eq. (jl]) by 1/2. 
Then the problem reduces to the Landau-Zener type 
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2G| 




FIG. 1. Energy spectra zbe^ for q = O.Tvr, O.Svr, O.OSvr from the outer curves to the inner. Tunneling 
happens only around t = r/2 and only for small q. 



Next, following Dziarmaga 16|, we shift the time t ^ t' 

t' 



rsing / t , 

— 1 - -(1 + cosg) 

1 + cos q \ T 



(6) 



Equation (jl]) is then transformed to 



d /u, 



'df 



-/i/ -1 \ /Uq 
-1 jjqt' j \Vq 



(7) 



where fig = 2(1 + cos q) / (r sin^ q). For a positive q, the transformed time t' runs in the range 



2 1 , 2 1 

-fi;' <t'< - — 



(8) 



smg " tang' 

For negative q, the signs of inequalities are reversed. Tunneling takes place during < 
t' < fi~^ because the off-diagonal elements of Eq. ([7]) are then larger than the diagonal 



elements in the absolute value. This latter time window is included in the above range of t' , 
Eq. dH]) for small |g| ^ 7r/3. We are therefore justified to extend the range of t' to (—00, 00) 
for small q because there is no transition in the extra time range then included. Then Eq. 
d?]) becomes equivalent to the Landau-Zener equation [19]. The Landau-Zener formula gives 
the transition probability Pq for each g, 



exp 



-27r( 



exp 



— vrrsm q 
2(1 + cosg) 



(9) 



For larger g, transitions do not happen. 

Let us next estimate the probability to obtain the ground state Pes in nonadiabatic 
quantum annealing. Similarly to the adiabatic case, we set the initial state as the ground 
state of the initial Hamiltonian, {uq,Vq) = (0, 1). Then we have 



GS 



n 

|g|<7r/3 





—TIT sin^ q 


^1 — exp 


2(1 + cosg)_ 



(10) 



It is convenient to take a logarithm, 

lnPGS= ' 

|g|<7r/3 





—TIT sin^ q 


^1 — exp 


2(1 + cos q)_ 



N 

— !■ — / dqln 



27r , 

iV 1 

27r ^ n 

n=l 



exp 



dqexp 



—vrrsm q 
2(1 + cosg) 

—riTiT sin^ g 
2(1 + cos g) 



(11) 



Since nonadiabatic transitions happen only for g ^ having small energy gaps, the following 
approximation is allowed, 

/ — nvrr sin^ g \ / nvrr 



Additionally, if we extend the integration range from (— tt/S, tt/S) to (—00, 00) in Eq. ffTTl) . 
the integration becomes Gaussian, which is valid for large r since the integrand is very small 
outside the range (— vr/S, tt/3). We then have a simple form of InPcs'- 



\nPns ^ —rJ24n= ^C(3/2), (13) 

TYWT W^l^ TxJt 

* n — 1 * 



n=l 



where, ({z) is the Riemann zeta function, 

oo ^ 

CW = E;ir' (14) 

n=l 

with C(3/2) ^ 2.61248. 

As is clear from the above equation, Pqs is a monotone decreasing function of the anneal- 
ing time r, which seems to be physically trivial. However, we repeat nonadiabatic processes 
many times, and obtain the ground state. Thus, the time necessary to obtain the ground 
state is the product of the annealing time r and the number of trials. On average, the 
number of trials converges to the inverse of the probability Pqs- For example, if Pes = 1/3, 
we have to repeat the process three times on average to reach the ground state. Therefore, 
we have to consider the quantity t/Pqs as the correct measure of computational complexity. 
We then again analyze, 

ln^ = lnr + ^C(3/2). (15) 

If we differentiate the above equation with respect to the annealing time r, we find that the 
above function takes its minimum at finite r, 

- = 3^(c(3/2)j . (16) 

We substitute this r into the right hand side of Eq. ( !T3|) and obtain the minimum value of 
the time necessary to obtain the ground state, 

mm^= ( — C(3/2)^ ^ 1.27732 x A^^. (17) 

It is also possible to estimate the standard deviation of the realization times, in addition 
to the average 1/Pgs, using Eqs. (fTTj) . ( IT6|) and (fT7|) as. 



a = v/e2(l + e-4/3)^(l + e^{l + e-^/^)^). (18) 

The time to reach the ground state may then be the product of r and l/Pcs + c""? which does 
not change the computational complexity (9(A^^). Therefore we have shown that nonadia- 
batic quatum annealing can succeed in the same time scale as the adtiabatic one, Eq. (jS]), 
for the one-dimensional transverse-field Ising model. 

In order to check the validity of a few approximations in the analyses around Eqs. ([9]) 
and (|T2l) . we numerically solved the Scrodinger equations for each mode (jl]) and calculated 



the time necessary for nonadiabatic quantum annealing t/Pqs for several single realization 
time r and the number of spins A^. Figure [2l shows that the total annealing time t/Pgs is 
propotional to A^^ for = 10, 20, ... , 100. This observation supports the result Eq. (1T7|) 
that the computational complexity is 0{N^). 
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1500 
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5000 10000 



FIG. 2. Numerical results for iV^ versus t/Pgs for = 10, 20, . . . , 100. 



III. SUMMARY AND DISCUSSION 

In the present paper, we proposed the method of nonadiabatic quantum annealing. We 
repeat trials of quantum annealing in the nonadiabatic time scale, and find the ground state 
among them. For the one-dimensional transverse-field Ising model, we analytically estimated 
the computational complexity of nonadiabatic quantum annealing and have shown that the 
computational complexity is the same as the adiabatic one (9(A^^). In order to the check 
the validity of our analyses, we numerically solved the Schrodinger equation, and found 
consistent results. 

As mentioned in Sec. [H this is a first step to open the door of nonadiabatic quantum 
annealing. We have to study other models, especially models having randomness j9l-lll|. It 
would also be useful to investigate the effects of thermal agitation in some details. 
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